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We calculate the entropy in a trapped, resonantly interacting Fermi gas as a function of temperature for a wide 
range of magnetic fields between the BCS and Bose-Einstein condensation endpoints. This provides a basis for 
the important technique of adiabatic sweep thermometry, and serves to characterize quantitatively the evolution 
and nature of the excitations of the gas. The results are then used to calibrate the temperature in several ground 
breaking experiments on 6 Li and 40 K. 
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The claims that superfluidity has been ob- 

served in fermionic atomic gases have generated great excite- 
ment. Varying a magnetic field, one effects a smooth evo- 
lution from BCS superfluidity to Bose-Einstein condensation 
(BEC) 0. In this Letter we use a BCS-BEC crossover 
theory to study the entropy S over the entire experimentally 
accessible crossover regime. Our goal is to help establish a 
methodology for obtaining the temperature T of a strongly in- 
teracting Fermi gas via adiabatic sweeps. This addresses an 
essential need of the experimental cold atom community by 
providing a temperature calibration for their ground breaking 
experiments |2j, |8j, |9J . In the process, we characterize quan- 
titatively the evolution of the excitations and show how their 
character evolves smoothly from fermionic to bosonic, and 
conversely. 

In adiabatic sweeps, the starting T at either a BEC or BCS 
endpoint is estimated from the "known" shape of the profile in 
the trapped cloud. Then, the temperature (near unitarity, say) 
is obtained by equating the entropy before the sweep to that 
in the strongly interacting regime after the sweep. Conven- 
tionally, the temperature scale which appears in the superfluid 
phase diagram |2, 8] involves an isentropic sweep between the 
unitary and the non-interacting Fermi gas regimes. The direc- 
tion of the sweep is irrelevant in these reversible processes. 
The important experimental phase diagrams plot the conden- 
sate fraction, N s /N near unitarity vs this Fermi gas-projected 
temperature, T e ff. 

In this paper, our thermodynamical calculations are used 
to relate the actual physical temperatures T to T e ff, where, in 
general, T is significantly greater than T P ff . A calculation of 
N s (T) is simultaneously undertaken fl(il.ll ill which provides 
an important self-consistency condition on the thermodynam- 
ics, since the same excitations appear in both, Moreover, a 
calculation of N s has to be done with proper attention paid to 
collective modes and gauge invariance 11211 . Here we address 
the various condensate fractions found experimentally 
(with emphasis on 6 Li), as a function of T e ff, in the experi- 
mental range of magnetic fields. 

Our work is based on the BCS-Leggett ground state 0,0] 
and its finite T extension 1 1 1 ]. Four different classes of exper- 
iments have been successfully addressed in this framework. 
These include (i) T w breathing modes experiments 
and theory II 3lll411 . (ii) radio frequency (RF) pairing gap ex- 



periments 1 9] and theory and (iii) T-dependent den- 

sity profiles 11711 . Finally, (iv) plots of the energy E vs T at 
unitarity 1 1 8 ] yield very good agreement with experiment and 
serve to calibrate the present thermometry. Two well-known 
weaknesses of the mean field approach (an underestimate of 
(3 and an overestimate of the inter-boson scattering length as 
in the deep BEC regime) should be noted. The first affects 
E{T) but not S(T). However, for the second we introduce a 
caveat: if the initial endpoint of sweep thermometry is suffi- 
ciently deep in the BEC regime (say, kpa < 0.3) the accuracy 
of the final temperature, we calculate for the unitary regime, 
could be improved by computing the initial S in the deep BEC 
regime using a pure-boson model with as set by hand to the 
Petrov result 1 19J . 

Because previous thermodynamic theories did not address 
unitarity, it has not been possible until now to arrive at a tem- 
perature scale in the experimentally interesting resonant su- 
perfluid regime. Carr et al. 1201 12111 calculated S at the BCS 
and weakly interacting, deep-BEC endpoints. The latter true 
Bose limit which they considered does not appear to be appro- 
priate to current collective mode experiments, O 01, which 
show that for physically accessible (i.e., near-BEC) 

fields, fermions are playing an important role. Thus, the BCS- 
Leggett ground state appears to be more appropriate than one 
deriving from Bose-liquid-based theory. Williams et al. I22I1 
calculated S for a BCS-BEC crossover theory using a mixture 
of noninteracting fermions and bosons 1 22] . This work, omits 
the important and self consistently determined fermionic ex- 
citation gap A which is an essential component for describing 
the thermodynamics of fermionic superfluids. 

Our thermodynamical calculations focus on this self- 
consistently determined A; they are based, for completeness, 
on a two-channel Hamiltonian II ll l23l 12411 . Here A appears 
in the fermionic dispersion = <J (e^ — [i) 2 + A 2 . (We 
define = h 2 k 2 /2m as the kinetic energy of free atoms, 
and the fermionic chemical potential.) Importantly, this A 
provides a measure of bosonic degrees of freedom. In the 
fermionic regime (\i > 0), A is just the energy required to dis- 
sociate the pairs and thereby excite fermions. At finite T, the 
closed-channel molecular bosons and the open-channel finite 
momentum Cooper pairs are strongly hybridized with each 
other, making up the "bosonic" excitations which contribute 
to thermodynamics. 



2 



Our many-body formalism has been described below the 
superfluid transition temperature r c |ll]. The parameter A, 
(when squared), is the analogue of the total number of par- 
ticles, in the simplest theory of BEC. Just as in BEC, there 
are two self-consistency conditions: (i) the effective chemi- 
cal potential of the pairs, fj, pa i r , is zero, for T < T c (as is 
that of the closed-channel molecular bosons f/, m b)> and, (ii) 
the number of pairs, reflected in A 2 (T) contains two addi- 
tive contributions representing condensed (A 2 C ) and noncon- 
densed (A 2 g ) pairs. The first condition implies that A(T) sat- 
isfies a BCS-like gap equation. Then, the condensate is de- 
duced, just as in BEC, by determining the difference between 
A 2 and A 2 g . In this approach the hybridized pairs have dis- 
persion flq — h 2 q 2 /2M*, with effective pair mass M*. 

We now extend this approach above T c . Our first equation 
represents the important defining condition on fi pa ir- that the 
inverse pair propagator (or T-matrix) t ~ 1 (Q)\q =0 — Z^pair, 
with (inverse) "residue" Z. While in the superfluid regions 
fipair and fi m f, vanish, in general, we have 

U efW + E 1 ~ 2 flf k) = Z ^r , CD 

where t/ e //(0) = U + g 2 / (2/i — v) involves the sum of the 
direct attraction U between open-channel fermions, as well as 
the virtual processes associated with the Feshbach resonance. 
Here f(x) is the Fermi distribution function. The determina- 
tion of the inter-channel coupling constant, g, and the mag- 
netic field detuning, v, is described elsewhere l25ll . as are the 
residues Z and Z II ill . The contribution from hybridized 
bosons will lead to a normal state excitation gap I9i. ll lll26l 12711 
or pseudogap (pg). This can be written in terms of the Bose 
distribution function b(x) as 

Kg = Z ' 1 zZ h ^ l i-^^- (2) 
q 

We use the local density approximation (LDA) throughout 
with a harmonic trap potential V(r). For notational simplicity, 
we omit writing V(r) in favor of /i(r) according to the LDA 
prescription: fj, — > fi(r) = /j, — V(r), where \i = ju(0). The 
total atomic number N = J d 3 r n(r), where 

n = 2n b0 + 2Z^ 1 ^b(n o - fi mb ) 
q 

+ 2j2[vl(l-f(E k ))+ulf(E k )} . (3) 

k 

Here n o — g 2 A 2 c /[(v — 2/j,(r))U} 2 is the density of con- 
densed closed-channel molecules, and u j, vj = [1 ± (ek — 
n(r))/Ei t ]/2. The total order parameter I23ll24ll26ll is given 
by A sc = A sc + \g\y/rib~o. 

To make progress, we numerically solve Eqs. Q-(|2) at each 
r for given /i and then self-consistently adjust \i via the total 
number constraint. Next we obtain the entropy S directly from 
the thermodynamical potential [28]. This potential contains 
fermionic contributions from bare fermions, ilf, and bosonic 




FIG. 1: Bosonic contribution to the thermodynamical potential. Here 
Go (G) and Do (D) are the "bare" ("full") propagators associated 
with the fermions and closed-channel molecular bosons, respec- 
tively, K and Q are four-momenta, and £/k,k' is the open-channel 
pairing interaction. 

contributions Vl . The latter is given by the sum of all possible 
ring diagrams shown in Fig.Q It can be easily shown that this 
fib is consistent with the self energy diagrams for the fermions 
and the molecular bosons. After regrouping, we see that S has 
two contributions, from fully dressed fermions (Sf) and from 
their bosonic counterpart {S ). The total entropy involves an 
integral over the trap, given by S = J d 3 r s(r) (and similarly 
for Sf and Sb), where 

S = Sf + Sb, 

s f = -2^[/ fc ln/ fe + (l-/ fe )ln(l-/ fe )], 

k 

s b = -E[ 5 9 ln6 9-( 1 + M ln ( 1 + M, ( 4 ) 

where fk = /(-Ek), and b q = 6(f2 q — [i m b)\ a relatively small 
contribution associated with the T dependence of f2 q has been 
dropped. The fermion contribution coincides formally with 
the standard BCS result for noninteracting quasiparticles [al- 
though here A(T C ) =^ 0]. And the bosonic contribution is 
given by the expression for non-directly-interacting bosons 
with dispersion £! q . These bosons are not free, however; be- 
cause of interactions with the fermions, their propagator con- 
tains important self-energy and mass renormalization effects. 

Figure [2] illustrates the behavior of S as a function of T 
obtained from our self-consistent equations, over the entire 
experimentally relevant crossover regime. The magnetic field 
is contained in the parameter 1/kpa, which increases with 
decreasing field. Here a is the s-wave fermionic scattering 
length, kp is the Fermi wavevector at the trap center, and 
ksTp = h 2 kp/2m is the noninteracting Fermi tempera- 
ture. Two important aspects of the fermionic contribution Sf 
should be noted. Generally, the fermions have a gap A in their 
excitation spectrum (which increases with decreasing field) 
and moreover this T dependent gap is inhomogeneous so that 
the fermions near the trap edge often behave as free particles 
at T > A. These quasi-"free" fermions change the T depen- 
dence of Sf from exponential to power law. They have also 
been seen in RF experiments |9, 15], as a free fermion peak in 
the spectra. 

We refer to Fig. [2] starting from the high field or BCS 
regime where S is linear in T. As the field is lowered to- 
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FIG. 2: (color online) Entropy per atom as a function of T for dif- 
ferent values of 1/fcpa from BCS to BEC in a harmonic trap. The 
dotted lines show an isentropic sweep between 1/kFa = 1 and uni- 
tarity. For comparison, we also plot S for an ideal Bose gas (dashed 
line). The 1/kFa = 3 curve lies below the dashed line at T > 0.8Tf 
reflects that M* ^ 2m. The inset plots the spatial profile of total en- 
tropy s (black curve) and its fermionic (s/, red) and bosonic (st, 
blue) component contributions at unitarity for T — T c /4. Here Rtf 
is the Thomas-Fermi radius, and T c = 0.27Tf. 



wards unitarity, Sf will vary as a low-T power law which 
is higher than linear. Simultaneously, the bosonic degrees 
of freedom emerge. Here one sees a T 3 / 2 power law from 
these excited bosons. At unitarity, bosonic effects dominate 
for T/T F < 0.05 or T/T c < 0.2. For an extended range of 
T < T c , the fermions and bosons combine to yield, S oc T 2 , 
which can be compared with the experimental power law 
T 2 - 73 . Finally in the near-BEC regime one sees an essentially 
pure bosonic T 3 / 2 power law in S at low T in the superfiuid 
phase. The relative contribution of the bosonic excitations, 
Sb/S evolves continuously from to 1 as 1/kpa increases 
from — oo to +oo. S becomes dominantly bosonic once fj, 
becomes negative. 

The bosonic T 3 / 2 power law found in the trap is the same 
as found for the homogeneous situation. Inhomogeneity effec- 
tively disappears here because the fermion-boson interactions 
lead to the self-consistent constraint that fj, pa i r — for the 
entire superfiuid region. This same disappearance of inhomo- 
geneity is found in Ref. 1 20] . This is different from a strictly 
non-interacting Bose system I22I1 (dashed line in Fig.|2j where 
one does not have a vanishing boson chemical potential below 
T c except at the trap center. The previous work of Ref. 1 20] is 
based on interacting but true bosons. The present situation is 
more complex since Cooper pair operators do not obey Bose 
commutation relations, (nor does the linear combination of 
Cooper pair and closed-channel boson operators). This sug- 
gests that a theory based on a true Bose liquid may not be 
appropriate for the fields that have been accessed experimen- 
tally. Moreover, if one were to contemplate contributions from 
the linearly dispersing Goldstone bosons, albeit within a more 



general ground state, their contribution, at unitarity, will not 
be as important as that from the edge fermions. 

To shed additional light on the component fermionic and 
bosonic contributions, in the inset to Fig.|2]we decompose the 
various terms in the entropy to reveal their spatial distribu- 
tions, for the unitary case at T = T c /4. It can be seen that 
the fermionic contribution s / (red curve) is limited to the trap 
edge, where A is small. By contrast, the bosonic contribution 
s;, (blue curve) is evenly distributed over the superfiuid region 
and rapidly decays at larger radii. 

Figure [2] provides a basis for thermometry in adiabatic 
sweep experiments. The vertical lines illustrate how to choose 
an initial temperature (TJ = 0.5Tf at point "A") with an initial 
value of 1 jkpa (=1) and use an isentropic sweep (represented 
by the horizontal line through "A" and "B") to obtain the final 
temperature (Tf = 0.28Tp at point "B") with the final value 
of 1/kpa (=0). It is most convenient to begin with either the 
BCS regime or BEC regimes, since here Tj can, in principle, 
be determined by fitting the density profiles. 

Figure [5] presents a plot of the superfiuid fraction fulfllll 
N s /N in the intermediate regime as a function of an effec- 
tive temperature T e g/Tp for different values of initial fields 
or 1/kpa. Here T e g is the temperature reached after an adi- 
abatic sweep to a BCS-like state. Based on experiment, we 
take the final state as l/k F a = -0.59 at 1025 G for 6 Li H 
and a noninteracting Fermi gas for 40 K 1 8]. 

The same vertical axis appears in the inset but with the 
physical temperature scale T, so that this figure provides a 
means of directly calibrating the temperature scale T e g which 
has been used in the important phase diagrams of 6 Li and 40 K. 
Figure [3] also provides a means of comparing the condensate 
fractions with those in the phase diagrams. For 6 Li at 900 G, 
with Teff/Tp = 0.2, 0. 1 and 0.05, the experimental condensate 
fractions are 0.0, 0.1 and 0.6. This should be compared with 
our calculated values, 0.006 0.36, and 0.73, respectively. For 
6 Li at 770 G, the condensate first appears at T eii /T F = 0.18, 
consistent with theory. From the values of kpa and T c at both 
ends (importantly, the latter can be read from the inset), one 
can easily see that the sweep from 770 G to 1025 G is still 
very far from a full BEC-BCS sweep. 

There are two reasons for the larger condensate fractions 
found theoretically for 6 Li. A calculation of the BEC-like 
density profiles shows that the noncondensed pairs inside the 
superfiuid region have a flat density distribution, which re- 
flects the vanishing of fi pa i r 1 17]. The superfiuid fraction 
extracted experimentally (assuming a Gaussian form for the 
noncondensed particles inside the condensate core) is, thus, 
underestimated, most notably around T c /2. In addition, ear- 
lier work 1 18] shows that when the system is treated as a non- 
interacting Fermi gas, T e ff will be underestimated whenever a 
condensate is present (at T < T c « 0.17T F at 1025 G). This 
suggests that theory and experiment can be brought into rather 
good agreement for the case of 6 Li. For 40 K, one has to ap- 
peal to non-adiabaticity and other complications of the sweep 
process to understand the small measured fractions. 

For this case, we emphasize temperature scales. For a full 
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FIG. 3: (color online) Superfluid density N a /N at different magnetic 
fields for 6 Li and 40 K as a function of the effective temperature, T e ff, 
measured in a near-BCS (at 1025G for 6 Li) or noninteracting Fermi 
gas (FG, 40 K) state accessed via reversible adiabatic sweeps of mag- 
netic field. The inset plots the same N s /N as a function of the phys- 
ical temperature at each field value. The system is not far from the 
resonance in these states. Also plotted in the inset is N S (T)/N at 
1025G and = 3 for 6 Li. The values of field and kpa were 

chosen based on Refs. 01 and H. For 6 Li, T F = 3.6/zK, so that 
l/k F a = -0.59, -0.26, 0.07, 0.38 for 1025 G, 900 G, 820 G and 
770 G, respectively. 

BCS to near-BEC [(£;Fa)finai — 1.7] adiabatic sweep 01 with 
initial T e g = Tj = 0.192f, the final reported temperatures in 
experiment |J and in theory are T f = 0.477> and 0.33T F , 
respectively. For the same (fcFa)finai but with T, = 0.17Tp, 
we find Tf w T c , in agreement with the observed sudden 
onset of a bimodal distribution in the density profile. Simi- 
larly for a sweep from a Fermi gas down to (fc^a)fi na i = 0.99 
with 7j = 0.061f, the experimentally quoted and theoreti- 
cally calculated Tf are 0.25Tp and 0.18Xi?, respectively. The 
experimental sweeps were not strictly adiabatic 1 8], so that the 
experimental Tf should serve as upper bounds. Our calcula- 
tions are more consistent with experiment, than if one had pre- 
sumed a T 3 power law for S in the BEC regime, from which 
one would infer Tf = 0.52Tp and 0.37Tp, respectively, ex- 
ceeding the upper bounds. 

Experimentally, 40 K gases J^l are prepared in the nonin- 
teracting limit, where, as a result of heating associated with 
an adiabatic sweep, low T is difficult to reach. By contrast 
6 Li gases |0, 01 are prepared in the BEC regime, so that 
higher condensate fractions of 80% and 95% have been re- 
ported |2l|4j,|9j] near unitarity via adiabatic cooling. Finally, 
we note that the rather large T c sa 0.1 7T F at 1025 G makes 
it hard to access the Fermi gas regime in 6 Li. This may be 
circumvented either by reduction in the size of N or Tp or, 
possibly, by sweeps to 528 G 1 29 ] . In 40 K, one avoids this 
problem altogether. 

Without knowing the temperature, measurements in this 



field cannot be directly compared to any theory. The present 
work presents a theory for the entropy S of a Fermi gas, at 
general accessible magnetic fields, which thereby calibrates 
T in various existing |2, 8, 9] and future experiments. 
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